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Abstract We investigated the environment of the infrared dust bubble N68 and searched 
for evidence of triggered star formation in its surroundings. We performed a multi wave- 
length study of the nebula with data taken from several large-scale surveys: GLIMPSE, 
MIPSGAL, IRAS, NVSS, GRS, and JCMT. We analyzed the spectral profile and the 
distribution of the molecular gas ( 13 CO J = 1 - and J = 3 - 2), and the dust in the en- 
vironment of the N68. The position-velocity diagram clearly shows that the N68 may be 
expanding outward. We used two three-color images of the mid-infrared emission to ex- 
plore the physical environment, and one color-color diagram to investigate the distribution 
of young stellar objects (YSOs). We found that the 24 /mi emission is surrounded by the 
8.0 jum emission. Morphologically, the 1.4 GHz continuum correlates strongly with the 
24 jim emission, and the 13 CO J = 1 - and J = 3 - 2 emissions correlate well with the 
8.0 yum emission. We investigated two compact cores located at the shell of the N68. The 
spectral intensity ratios of 13 CO J = 3- 2toJ=l-0 range from 5 to 0.3. In addition, 
young star objects, masers, IRAS, and UC HII regions distribute at the shell of bubble. 
The active region may be triggered by the expanding of the bubble N68. 



Key words: infrared: stars — stars: formation — ISM: bubbles — HII regions 



1 INTRODUCTION 

There are many signatures of star formation, for example, outflow/inflow, dark cloud and HII re- 
gion, which can be used to investigate the process of star formation. Also there are m any different 
kinds of interactions which can trigger the st ar formation, such as cloud-cloud col lision dLi & Wangl 

" l2012h ar 



20121). supernova explo s ion (|Xu et all 1201 l a b), bubble expansion dZhang & Wangl I20T2I) and so on. 
Church well et al. I (I2006L 1200717 detected and cataloged about 600 mid-infrared dust (MIR) bubbles be- 



tween longitudes -60° and +60°. The bubbles have bright 8.0 jim shells that enclose bright 24 jim 
interiors. The infrared (IR) dust bubbles may be produced by exciting O- and/or B-type stars, which is 
located inside bubble. The ultraviolet ( UV) radiation from exciting stars may heat dust and ionize the 
gas to form an expanding bubble shell dWatson et alll2008b . which is known as "collect-and-collapse" 
process. This process can trigger a massive star formation near the shell clumps. Therefore, the bubbles 
present an important opportunity to study the interaction between the HII regions and molecular clouds. 

A few individual bubbles have been studied well, such as N49 and S51. Ther e are many models and 
observ ations to explain the dusty wind-blown bubbles, such as bubble N49 of lEverett & Churchwelll 
(I2010h . Recently, we reported an expan ding bubble S5 1 shown with a shell and a front side, e mploying 
13 CO and C 18 J = 1 - emission lines dZhang & WangLl2012h . lBeaumont & Williams! (l2010l) reported 
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CO J = 3 - 2 maps of 43 Spitzer identified bubbles. IWatson et al.l d2008[) present an analysis of wind- 
blown, parsec-sized, mid-infrared bubbles and associated star formation. 

To add the investigate d bubble example s , we se lected the IR dust bubble [CPA2006] N68 (hereafter 
N68) from the catalog of IChurchwell et alJ (120061) . N68 is a complete or closed ring and centered on 
Galactic /=35.654, £=-0.062 (or R.A.(J200 0) = 18^56^5^.70 DEC (J2000) = 02°26 , 01 ,, .0). It has a 
distance of 10.6 kpc, which was obtained by I Anderson & Banial ([2009) based on comparing the velocity 
of the ionized gas with the maximum velocity of HI absorption, and one looking for HI absorption at the 
velocity of molecular emission. In addition, its size is 34 pc x 17 pc and the eccentricity of the ellipse 
is 0.72. 

In this work, we mainly present a multiwavelength study of the environment surrounding the IR 
dust bubble N68. We explore its surrounding interstellar medium (ISM) and search for signatures of 
star formation. The observations and data are described in Sectj2l the results and discussions about the 
bubble N68 environment are presented in Sect|3j Sect|4] summarizes the results. 

2 DATA 



Table 1 YSOs around the bubble N68 



Source 


Name 


R.A.(J2000) 

h m s 


DEC.(J2000) 


IRAS 


IRAS 18538+0222 


18 56 21.30 


+02 26 37.00 


IRAS 


IRAS18538+0216 


18 56 23.50 


+02 20 38.00 


IRAS 


IRAS 18540+0220 


18 56 35.60 


+02 24 54.00 


CHII 


G35.590-0.025 


18 56 22.67 


+02 21 14.58 


EGO 


G35.040-0.470 


18 57 03.30 


+02 21 50.00 


Maser 


G35.578-0.031 


18 56 22.55 


+02 20 28.10 


UCHII 


G35.578-0.031 


18 56 22.64 


+02 20 26.30 



We analyze d IR and millimeter wavelength data extracted from several l arge-scale surveys includ- 
ing GLIMPSE dBeniamin e t all 120031: IChurchwell et all l2009h. MIPSGAL dCarev et alll2009h . IRAS 



dNeugebaueret allll984l) . NVSS (ICondon et allll998h . GRSf 1 ! and JCMT0. 



GLIMPSE is a MIR survey of the inner Galaxy performed with the Spitzer Space Telescope. We 
used the mosaicked images from the GLIMPSE and the GLIMPSE Point-Source Catalog (GPSC) in 
the Spitzer-IRAC (3.6, 4.5, 5.8, and 8.0 /mi). IRAC has an angular resolution of between 1.5" and 1.9" 
dFazio et all 120041: IWerner et all 12004 . MIPSGAL is a survey of the same region as the GLIMPSE, 
using the MIPS instrument on Spitzer. The MIPSGAL resolution is 6" at 24 /mi. The 60 and 100 
yum fluxes of the IR Astronomical Satellite (IRAS) were employed to estimate the Lyman continuum 
ionizing flux. The NRAO VLA Sky Survey (NVSS) i s a 1.4 GHz continuum survey covering the entire 
sky north of -40 deg declination dCondon et allll998|) . 

The 14-m GRS has a full width at half maximum (FWHM) beam size of ~ 46" at -110 GHz for 
13 CO J = 1 - transitions. Each field of the GRS in the 13 CO J = 1 - line comprises spectra on a fully 
sampled 22" grid. The intensities are on a T* A antenna temperature scale. To convert this to main beam 
temperature, is divided by the m ain beam effic i ency o f 0.48. The velocity resolution of the data is 
0.25 km s _1 (0.22 km s _1 sampling) dJackson et al.l.l2006h . The 15-m JCMT has a FWHM beam size of 
~ 14" at ~330 GHz for 13 CO J = 3 - 2 transitions. The correlator was configured with 4096 channels 
over a 250-MHz bandwidth, which provided a velocity resolution of -0.055 km s" 1 per channel. The 



1 This publication makes use of molecular line data from the Boston University-FCRAO Galactic Ring Survey (GRS). The 
GRS is a joint project of Boston University and Five College Radio Astronomy Observatory, funded by the National Science 
Foundation under grants AST-9800334, AST-0098562, AST-0100793, AST-0228993, & AST-0507657. 

2 The James Clerk Maxwell Telescope is operated by the Joint Astronomy Centre on behalf of the Science and Technology 
Facilities Council of the United Kingdom, the Netherlands Organisation for Scientific Research, and the National Research 
Council of Canada. 
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Fig. 1 Mid-IR emission of the IR dust bubble N68. Left: Spitzer-IRAC three-color image (3.6 
yum = red, 4.5 yum = blue, and 8.0 yum = green); Right: Spitzer-IRAC and Spitzer-MIPSGAL 
three-color image (4.5 yum = blue, 8.0 fim = green, and 24 fim = red). 



spectra of 13 CO J = 1 - and J = 3 - 2 have been smoothed to a velocity resolution of 0.50 and 0.44 km 
s" 1 , respectively. 

The 13 CO J = 1 - and the 13 CO J = 3 - 2 data cubes were processed with CLASS and GREG in 
the GILDAS software package^. 

3 RESULTS AND DISCUSSIONS 

3.1 The IR dust distribution of the N68 

Two Spitzer three-color images of the N68 were shown in FigQ] Both figures clearly illustrate the 
photo-dissociation region (PDR) visible in 8.0 fim (in green) emission, which originates mainly in the 
poly cyclic aromatic hydrocarbons (PAHs). Due to these large molecules inside the ionized region are 
destroyed, the PAH emission delineates the HII region boundar ies, and molecular clouds inside are 
excite d in the PDR by the radiation leaking from the HII region (IPetriella et all l2010l : IPomares et all 
120091) . The 24 yum emission (right in FigQ]), appearing just in the northwest part of the bubble N68, 
corresponds to hot dust. It is likely that O- and/or early B-type stars produced the bubble shell of this 
HII region, with hot dust located inside the bubble. The 3.6 yum emission (in red, left panel in FigQ]) 
shows the positions of different stars. 

In addition, the molecular gas of the bubble shell exhibits several clumps and filaments along the 
PDR shell seen from FigQ] The distribution and morphology of this material suggests that a collect- 
and-collapse process may be occurring. 

3.2 The continuum emission and clump structure 

The 1 .4 GHz continuum NVSS is full of the entire bubble N68 (left panel in Figj2]), particularly saturated 
at the northwest part. The area, saturated for the NVSS, morphologically correlates with the 24 yum 
emission in FigQ] We suggest that exciting stars, which are radiating strong UV radiation, may exist in 
this region. Possibly because the obstruction from molecular cloud, the hot dust develops and permeates 
into the southeast part of N68. 



3 http://iram.fr/IRAMFR/GILDAS/ 
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Fig. 2 Integrated velocity contours of 1.4 GHz continuum (Left panel) and 13 CO J = 3 - 2 
(Right panel), superimposed on the GLIMPSE 8.0 fim color images. The contour levels are 
from 0.16 to 1.45 by 0.16mJy/beamfor 1.4 GHz continuum, while from 6.08, 12.17 to 121.70 
by 12.17 K km s" 1 for 13 CO J = 3 - 2. The integrated velocity range is from 38.5 to 68.8 km 
s -i f or 13£Q j = 3 _ 2. The symbols "a", and indicate the positions of the maser 
(H2O and OH), the IRAS point source, and the EGO source, respectively. The beam size is 
shown in black filled circle. Note that the 1.4 GHz continuum are saturated with oblique line 
in the image. 



The 13 CO J = 3 - 2 contours (Figj2j) from JCMT data show several clumps, which may be the birth 
place of YSOs. These contours are almost distributed on the 8.0 /zm shell of the bubble. At the southwest 
of the N68, there are two strong cores which will be investigated in the followed Sect l3.4[ 

The 1.4 GHz continuum flux is about 1.27 Jy within this bu bble. The number o f stellar Lyman 
photon, absorbed by the gas in the HII region, follows the relation dMezger et all Il974l) 

[^] = 4.761 x 10 48 • a(v, T e )~ y ■ [-^-] 01 • f^]" 045 • [^] • [^-f, (1) 
s~ l GHz K Jy kpc 

where a(v, T e ) is a slowly varying function tabulated by lMezger & Henderson! d 1967k : for effective tem- 
perature of exciting star T e ~ 33000 K and at radio wavelengths, a(v, T e ) ~ 1. Finally, we obtain Lyman 
continuum ionizing photons flu x logNr ~ 48.7 9 from nebula. Assuming the exciting stars belong to 
09.5 star with logA^ = ~ 47.84 dPanagialll973h . we suggest there should be about 8.89 exciting stars to 
ionize the bubble. And there are presumably many other non-IR excess sources within the bubble that 
could also be O- and/or early B-type stars. 

3.3 The distributions of the 13 CO J = 1 - molecular clouds 

3.3.1 Opacity, column density, and mass 

We used 13 CO J = 1 - to derive the opacity, column density, and mass of nebula within the box shown 
in Figj3] The size of the box is 13.3'xl3.3'. The velocity range is from -38.5 to -68.8 km s" 1 . The box 
is centered on R.A.(J2000) = 18^56 m 25 5 .70, DEC.(J2000) = 02°26 , 01 ,, .0. 

Here we assume that the nebula is in local the rmodynamic equilibrium ( LTE). We used the t heory 
of the radiation transfer and molecular excitation dWinnewisser et all 1 19791: iGarden et all Il99lb . The 
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Fig. 3 Left: integrated velocity contours of the 13 CO J = 1 - emission superimposed on the 
GLIMPSE 8.0 yum color image. The contour levels are from 14.23 to 61.66 by 4.74 K km 
s" 1 , and the integrated velocity range is from 38.5 to 68.8 km s" 1 . The dotted ellipse shows 
the PAH morphology of bubble. These numbers give the chosen positions of seven sources. 
The YSOs candidates of class I, II are indicated with symbols "x", and "a" respectively. The 
symbols "a", and indicate the positions of the maser (H 2 and OH), the IRAS 
point source, and the EGO source, respectively. The beam size is shown in black filled circle. 
Right: the position- velocity diagram of the 13 CO J = 1 - emission along the dashed line in 
left panel. The contour levels are from 0.35 to 3.18 by 0.36 K. 



i * * * * i * * * * i 




[5.8]-[8.0] 



Fig. 4 The color-color diagram between the [5.8-8.0] and [3.6-4.5] of the GLIMPSE. Class I, 
II and III stars of different evolutionary stage have been indicated in the diagram. 
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Fig. 5 13 CO J = 1 - (red line) and J = 3 - 2 (gray line) spectra of sources 1 to 7, whose 
positions are represented in Fig [3] with number 1 to 7. All indicates that the two spectra is the 
spectral average of the bubble N68. 
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CO J = 1 - opacity is gained from the equation: 



r( 13 CO) = -ln[l - TMB ^ C0 \ [eX p( 5 ^ - I)]" 1 - 0.16}- 1 ], 

lex 



(2) 



where the T ex is excitation temperature. The 13 CO J = 1 - column density is obtained from the equation: 

T ex fr( l3 CO)dv 



N( l3 CO) = 2.4 x 10 



14 



l-exp(-5.29/T ex )' 



(3) 



We assume that the 13 CO J = 1 - is optically thin, t he exc itation temperature is T ex - 20 K and the 
[H2/ 13 CO] abundance ratio is 7 x 10 5 dFrerking et all fl982b . The molecular hydrogen column density 
N(H2) was then calculated. The molecular cloud mass is estimated at M#68 ~ 4.6 x 10 5 M© including 
the entire bubble N68 shell at a distance of D N ^ ~ 10.6 kpc. 



3.3.2 Clump structure 

Fig (3] shows integrated velocity contours of the 13 CO J = 1 - emission superimposed on the GLIMPSE 
8.0 yum color image. The integrated velocity range is from 38.5 to 68.8 km s" 1 . Generally, there should be 
similar sketch between th e distributions of the millimeter emission and Mid-IR emission dPetriella et all 
l2010l:lzhang et all 1201 2h . Comparing the morphologic distribution (Fig (3]), there are many correlations 
between the contours and the 8.0 pm emission, however, their peak positions deviate slightly from each 
other. In addition, the nebula is clumped and extended outside the bubble. Therefore we have selected 
seven sources (1 ~ 7; indicated in FigO to investigate their spectral information in followed Sect l3.3.4l 
Along the dash-straight line in Fig0 we show the position- velocity diagram. Fig(3] shows that the 
components of sources A, B, C, D and E are correlative with the N68, while the component of source 
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Fig. 6 Integrated velocity contours of the 13 CO J = 1 - and J = 3 - 2 emission superimposed 
on the GLIMPSE 8.0 and 24 jim color images, respectively. The contour levels are from 14.23 
to 61.66 by 4.74 K km s" 1 for 13 CO J = 1 - 0, while from 6.65, 13.30 to 119.68 by 13.29 K 
km s" 1 for 13 CO J = 3 - 2., and the integrated velocity range is from 38.5 to 68.8 km s" 1 . The 
dotted ellipse shows the PAH morphology of bubble. These numbers give the chosen positions 
of three sources. The YSOs candidates of class I, II are indicated with symbols "x" and "a". 
The symbols "a", unfilled star, and indicate the positions of the maser (H2O and 
OH), the IRAS point source, the compact HII region, and the UC HII region, respectively. 
The white circle indicates the position error of the IRAS point source. 



F is not. The components of sources A, B, C, D and E are made up of a morphology of letter "N" in 
the position-velocity diagram (right panel in Fig©. Source C inside bubble N68 has only one velocity 
component. And from source B (in the shell of the bubble) to source D (in the shell of the bubble), the 
position- velocity diagram dose not show the l etter "U" shape. This is not followed by the model of an 
expanding bubble inside a turbulent medium dArce et allEoTTh . However, considering the eccentricity 
of the bubble and positions of sources A, B, C, D and E, we also suggest that the bubble N68 may be 
expanding outside. The systematic velocity of the N68 is -58 km s" 1 (the source C inside the bubble), 
and the velocities are ~59 and -57 km s" 1 for sources A and B outside the bubble, respectively. Also, 
the velocities are ~63 and -53 km s" 1 for sources B and D in the shell of the bubble. So the bubble N68 
may be expanding with a velocity of ~5 km s" 1 along the line of sight relative to source C. 

3.3.3 The young stars distribution 

The color-color diagram (Figj4j) shows the distribution of class I, II and III st ars. Here we only consider 
these sources with detection in four Spitzer-IRAC bands dHora et all l2008|) . In addition, to look for a 
relation between the young stars and the bubble formation, we draw the class I and II stars on the 8.0 
yum emission image of the N68 (FigO. We found that class I and II stars almost are distributed on the 
shell of the bubble. The distribution of class I and II stars is also correlative with 13 CO cloud. We also 
report the distributions of the IRAS, the extended green object (EGO), the maser and the UC HII region. 
These star formation tracers are located near and on the shell of nebula. It is possible that the shell of 
the bubble is the birth place of young stars, which may be the result of interacting effect between the 
HII region and molecular cloud, and be the production of a collect- and-collapse process. 
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Fig. 7 13 CO J = 1 - (red line) and J = 3 - 2 (gray line) spectra of compact cores 
M (R.A.(J2000) = 18^56 m 33^73, DEC.(J2000) = 2°20 , 49 ,, .81) and N (R.A.(J2000) = 
18*56 m 23 5 .07, DEC.(J2000) = 2°20 , 33 ,, .44), whose positions are also represented in FigEl 



3.3.4 The spectra of the 13 CO J = 1 - and J = 3 - 2 

FigH shows the 13 CO J = 1 - and J = 3 - 2 spectra of sources 1, 2, 3, 4, 5, 6, and 7 on the shell of 
the nebula. The average of all spectra of the whole bubble is exhibited (indicated with All) in panel of 
Fig |3 This shows that the spectra of the 13 CO J = 1 - has a wider width and stronger intensity than 
those of the J = 3 - 2. At the sources 2 and 5, however, the spectra of the 13 CO J = 3 - 2 has a stronger 
intensity than that of the J = 1 - 0. The 13 CO J = 1 - spectrum of source 4 does not have the -61 km 
s" 1 component relative to the J = 3 - 2. In addition, from the contours, the 13 CO J = 3 - 2 distribution 
(FigO is clumpier than that of the J = 1 - (Fig|2|). Therefore, the difference between the 13 CO J = 1 - 
and J = 3 - 2 is very much, and needs further investigation. 

The intensity ratio of the 13 CO J = 3 - 2 to J = 1 - is Rl\ n n = 0.46, 3.53, 0.38, .70, 2.69, and 0.46 
for sou rces A//, 2, 3, 4, 5, and 6, respectively. Our results are consistent with that of iMinamidani et al.l 
d201ll) . whose R l 3 3 _ 2/l _ ranges from 3.6+0.9 to 0.24+0.07. This is suggestive that 13 CO J = 1 - and J 
= 3-2 transitions have different critical densities, so the ^3 3 2/1 _ traces the excitation of molecular gas. 

3.4 Two compact cores in the shell of the N68 

There are two compact cores M (R.A.(J2000) = 18 /l 56 m 33'.73, DEC.(J2000) = 2°20 , 49 ,, .81) and N 
(R.A.(J2000) = 18 /z 56 m 23^.07, DEC.(J2000) = 2°20 , 33 ,, .44), located at the southwest shell of N68. 
Fig|6l shows the integration contours of the 13 CO J = 1 - and J = 3 - 2 emission. The cores M and N 
are located at the peaks of integration contours of the 13 CO J = 3 - 2 emission. The peak of the 13 CO J = 
1-0 contours is located at core N, but it does not appear for core M. In addition, the difference between 
the 8.0 and 24 jim emission is obvious. The present of the 24 jim emission is suggestive that the hot dust 
has plunged into the bubble shell. 

On the other hand, there are many star formation tracers located near core N. We marked the posi- 
tions of a H 2 maser, a OH maser, a IRAS point source, a compact HII region and a UC HII region in 
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Fig 16) The H2O maser, the OH maser, the IRAS point source and the UC HII region are located at the 
core N. There are many velocity components for the H2O maser and the OH maser. We just list three of 
their strongest components. The fluxes are 20.15 Jy at 48.90 km s" 1 , 7.25 Jy at 49.45 km s" 1 and 2.80 Jy 
48.35 km s _1 for H 2 maser emission, whil e 9.90 Jy at 50.00 km s" 1 4.40 Jy 46.05 km s _1 and 3.70 Jy 
at 47.36 km s" 1 for the OH maser emission dForster & CaswelllJl999h . These are suggestive that core N 
is a strongly active massive star forming region. 

We also exhibit the spectra of cores M and N in Fig(7J The difference of spectral profiles between 
13 CO J = 1 - and J = 3 - 2 emission is obvious, the 13 CO J = 3 - 2 has double profiles for cores M 
and N, but the 13 CO J = 1 - lacks of the blue profile. For core N, the 13 CO J = 1 - has a main beam 
temperature of ~7 K at -53 km s" 1 , yet the 13 CO J = 3 - 2 shows a flat spectrum of absorption. This 
discrepancy is needed to further investigate. 

4 CONCLUSIONS 

We have investigated the environment of the IR dust bubble N68 with the 13 CO J = 1 - and J = 3 - 2 
emission lines, the 1.4 GHz continuum, and the IR bands emission. The main results can be summarized 
as follows. 

The morphology of the molecular line emission ( 13 CO J = 1 - and J = 3 - 2) and the associated 
velocity signatures are consistent with the shell structure seen from the 8.0 fim images. The expanding 
shells in N68 are suggestive of triggered massive star formation. There is evidence for dense clumps 
coincident with the shells. The expanding shells have an expansion speed of ~5 km s" 1 . However, it is 
not clear whether HII region is the driving mechanism of the shell expansion. 

We have presented a study about compact cores M and N located at the shell of nebula. Compact 
core N is associative with the H 2 maser, the OH maser, the IRAS point source, and the UC HII region. 
The spectral velocity components of the H 2 maser and the OH maser indicate that the star formation 
region is strongly activate. However, the deep-going properties of the compact cores M and N are needed 
to further investigate. 

We also used the GLIMPSE, and the MIPSGAL survey data to analyze the YSOs and the warm dust 
distribution around the bubble N68. We identified class I and II YSOs using the [5. 8] -[8.0] versus [3.6]- 
[4.5] relation, and correlated their distribution relative to the PDR, which we assume to be associated 
with and surrounding an HII region. We find that the N68 appears to have a significant number of YSOs 
associated with their PDRs, implying that triggered star formation mechanisms acting on the boundary 
of the expanding HII region. 
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